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Abstract 

In this comment we question some arguments presented in astro-ph/0512048 
to refuse the presence of an singular mass surface layer. In particular, 
incorrect expressions are used for the disk’s surface mass density. We 
also point out that the procedure of removing the descontinuity on the 
z = 0 plane with a region of continuous density gradient generates 
other two regions of descontinuities with singular mass surface layers 
making the model unrealistic. 


1 Introduction 

In [Ij the authors make some comments about their previous work ;2] and 
also reply to certain issues that have been raised by some authors El El El 
concerning the nature of the matter distribution and the asymptotic be¬ 
haviour of the spacetime in their model. In particular, .3] pointed out that 
the use of a \z\ function in the solution of the field equations introduces an 
additional disk on the 2 = 0 plane and later 0 showed that the disk was 
made of exotic matter. Consequently, [T] examine the question in several 
ways and present some arguments to dismiss the existence of a physical mass 
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layer present on the z = 0 plane. We would like to notice some inconsis¬ 
tencies in their paper. The first is related to the expressions used for the 
surface density of the singular disk, and is presented in Sec. |21 In Sec. 01 
we question the smoothing procedure discussed in Section 4 of IJ. We find 
that one cannot eliminate disk singularities, moreover one introduces new 
singularities that make the model unrealistic. 


2 The disk’s surface mass density 


In [I] we calculate the distributional energy-momentum tensor Q% due to the 
introduction of an absolute value of z in the metric functions. The resulting 
non-zero components read [5. 
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where all quantities are evaluated on z = 0. Since Q% is non-diagonal, in 
order to have obtain the physical variables of the disk, we need to put the 
energy-momentum tensor in its canonical form. To do that we solve the 
eigenvalue problem: Q^ b = A£“. We thus find 


A+ = — ± 


T , VD 


T = Q\ + QZ, D = (Q\ - + 4 Q^Qf, 
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and using Eq. (fTal - Jldl) result in 
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As the discriminant is always negative, the eigenvalues and corresponding 
eigenfunctions are complex conjugate. If V a and W a denote the timelike 
and spacelike real eigenvectors, the canonical form of the energy-momentum 
tensor is 

Qab = <rV a V b + p v W a W b + K{V a W b + W a V b ), (5) 
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where a = T/2 is the surface density, p v = — T /2 denote the azimuthal 
stresses and k = y/—D /2 is the heat flow in the azimuthal direction |B]. Thus 
the expression for the surface density, to order G 1 , is given by a = — v >z , or, 
using the relation i/ z = —N iT N iZ /r we get 


a = 


N r N z 

r 


( 6 ) 


On the other hand, [T take following expressions for the surface energy 
density (Equations (15) and (16) of the paper) 


NN Z NN Z N r N z 

a = - 2 - V * = - 2 + -’ ( 7 ) 

'j* Z. fyt ry* 

which they integrate over the surface and compare with the volume integral 
of their continuous mass density distributions. But Eq. 0 is only the Q\ 
component Eq. (HaJ) to order G 1 . Due to the non-diagonal form of the 
energy-momentum tensor, the Q\ component solely does not determine the 
surface density, but there is also a contribution from the Q^ component Eq. 
m even to order G 1 . Thus, the correct expression for cr that should be 
used is Eq. The integration of this equation over the surface would 
result in half of the value for the mass derived from the volume integral of 
the continuous mass distribution, since the authors themselves comment on 
footnote 7 that the two terms in Eq. 0 contribute equally. 

It is important to stress that the physical variables of the mass layer on 
the z = 0 plane are not directly given by the principal diagonal terms of 
the distributional energy-momentum tensor, since the non-diagonal terms 
are non-negligible. This is further an example of how the non-linearity of 
General Relativity can manifest even at Newtonian level. 


3 The smoothing procedure 

Another approach used in [Tj to examine the presence of a singular mass 
layer on the symmetry plane was to smooth the solution over an interval 
that includes the z = 0 plane. This was achieved by the choice of cosh(fc ri z) 
functions to span the symmetry plane in the interval — zo < z < zq. For 
|z| > zo, the original solution with exponentials was used. This requires 
that the functions N and also A( z match at \z\ = zq. If this last condition 
is not satisfied new matter is added. We note that this is also true for the 
Newtonian gravitational potential, discontinuity of the first derivatives “add 
matter” whose density can be computed via Poisson’s equation. We argue 
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Figure 1: Impossibility of matching the functions Eq. (I8a,l) m and their 
derivatives, (a) Continuity of N. (b) Continuity of N z . 


that the above mentioned matching cannot be done without adding new 
matter. 

Let us take the function N as follows: 


Ni = - Y, C ln k n e knZ r Ji(k n r), z < -z 0 , (8a) 

n 

N 2 = - ^2 Cinkn cosh(k n z)rJi{k n r), -z 0 < z < z 0 , (8b) 

n 

N 3 = -^2C 3n k n e~ knZ rJi(k n r), z > z 0 , (8c) 

n 


and impose the conditions N\{—zq) = N 2 (—zo) and IVi )Z (— zq) = 
Using Eq. m-m we obtain 


^ k n rJ\ (k n r) C 2n cosh(fc n zo) - C ln e knZ ° 

n 

^2 kn r Ji{k n r) ~C 2n sinh(fc n 2 :o) - C ln e~ knZ ° 


= 0 , 

= 0. 


N 2 ,z(~ z o)- 


(9) 

( 10 ) 


Since the Bessel functions are linearly independent, the terms in brack¬ 
ets must vanish identically. Subtracting Eq. © from Eq. enj results in 
C 2n e knZ ° = 0, which is only satisfied in the real domain if all C 2n = 0. Thus 
it is not possible to match N and N tZ simultaneously. If we demand con¬ 
tinuity of N\, N 2 and N 3 at \z\ = zq we obtain the following conditions: 


Cin — C 3n , and C 2 n 


Cm 


o knZ o 


cosh(fe n zo) 


( 11 ) 
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On the other hand, demanding continuity of 
results in 

Ci n — flirt ■ and C 2 n — Cm 


,z, N 2 , z and N 3:Z at \z\ = z 0 


g knZ 0 

sinh(A; n zo)' 


( 12 ) 


Fig. Ufa) sketches what happens if we impose continuity of the functions 
and Fig. mb) if we impose continuity of the derivatives. The only way the 
functions Eq. (j5al) -(jSc)) and their derivatives could be matched is if an extra 
set of constants were inserted into Eq. (|8b|) 


N2 = -Y: k n(C 2 n cosh (k n z) + Cf* tra )r J x (k n r ), (13) 

n 

but then this would not be a solution of 

Nrr + N ZZ -^L = 0. (14) 

r 

Assuming conditions Eq. m hold, the jump of the derivatives of Eq. 
(I8al)-(l8cl) with respect to z evaluated at \z\ = zq are given by 

N 2 , z (-zo) ~ Nx^-zo) = N 3tZ (z 0 ) - N 2tZ (zo) = ^ Clnfc ^ l(fc " r) - ( 15 ) 

cosn.[K n zo j 


The same kind of descontinuities in the derivatives also appear when a \z\ is 
introduced in the solution. Thus they introduce additional layers of matter 
now located on the z = ±00 planes that makes the solution unrealistic. 
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